Derealization and spin-wave dynamics in ferromagnetic chains with long-range 

correlated random exchange 
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We study the one-dimensional quantum Heisenberg ferromagnet with exchange couplings ex- 
hibiting long-range correlated disorder with power spectrum proportional to l/k a , where k is the 
wave- vector of the modulations on the random coupling landscape. By using renormalization group, 
integration of the equations of motion and exact diagonalization, we compute the spin-wave localiza- 
tion length and the mean-square displacement of the wave-packet. We find that, associated with the 
emergence of extended spin- waves in the low-energy region for a > 1, the wave-packet mean-square 
displacement changes from a long-time super-diffusive behavior for a < 1 to a long-time ballistic 
behavior for a > 1. At the vicinity of a = 1, the mobility edge separating the extended and localized 
phases is shown to scale with the degree of correlation as E c oc (a — l) 1 ^ 3 . 

PACS numbers: 75.10.Jm, 75.30.Ds, 75.50.Lk 



I. INTRODUCTION 



It is well established that one-electron eigenstates in 
chains with uncorrelated disorder are exponentially local- 
ized Q. However, several theoretical investigations have 
shown that a series of one-dimensional versions of the 
Anderson Model exhibit a breakdown of the Anderson's 
localization induced by internal short-range correlations 
on the disorder distribution, including hopping and on- 
site energy correlations Q or just on-site energy correla- 
tions H Q] . Recently, it has been demonstrated that the 
one-dimensional Anderson model with long-range corre- 
lated disorder can display a phase of extended electronic 
states H [?). These results have been confirmed by mi- 
crowave transmission spectra of single-mode waveguides 
with inserted correlated scatterers j7j. The influence of 
long-range disorder on the electron motion in two di- 
mensions has been recently investigated [|| and a two- 
dimensional layered Anderson model with long-range cor- 
related disorder has been shown to exhibit a Kosterlitz- 
Thouless-type metal-insulator transition [B|. 

It is also known that the magnon equation of mo- 
tion for ferromagnetic spin chains with uncorrelated ran- 
dom nearest-neighbor exchange couplings can be exactly 
mapped onto an electronic chain with a particular form 
of off-diagonal disorder where the random hopping in- 
tegrals appear correlated in pairs [|l0[ [ll]]. By using a 
perturbation approach combined with a scaling hypoth- 
esis it was demonstrated JTo| that the singularities of the 
density of states and the localization length of a random 
ferromagnetic Heisenberg chain depend on the distribu- 
tion of exchange couplings. For uncorrelated random ex- 
change couplings J G [0, 1] with a probability distribu- 
tion function P(J) = (1 - 5)J~ S , it was demonstrated 
that for S < — 1, which implies < l/J > and < 1/J 2 > 



finite, the density of states p(E) diverges as 1/E 1 / 2 and 
the localization length as 1/E. Distinct regimes emerge 
when either of the above two moments diverges, thus 
generalizing other studies E3. Moreover, in order to 
elucidate doubts raised in the literature [ fl3| , the spin- 
wave dynamics of the one-dimensional Heisenberg fer- 
romagnet was investigated ]IT| for the above power-law 
distribution of exchange couplings J. The mean-square 
displacement of the wavepacket asymptotically displays 
super-diffusion dynamics [cr(t) 2 oc t 3 / 2 ] for weak disorder 
(S < 0), diffusive behavior \a(t) 2 oc t 1 } for S = 0, and lo- 
calization [<r(t) 2 oc constant] for strong disorder (5 > 0). 
Therefore, the uncorrelated Heisenberg ferromagnet with 
random exchange couplings can display a super-diffusive 
dynamics if < 1/J > is finite. In all three cases, a bal- 
listic regime [u(t) 2 oc t 2 } emerges for initial times. By 
using the transfer matrices technique, the singularities 
of the density of states and localization length were veri- 
fied. The super-diffusive behavior is closely related to the 
one found in the random-dimer version of the Anderson 
model |J. 

In this paper, we study the nature of the spin-wave 
modes of a quantum Heisenberg ferromagnetic chain with 
long-range correlated random exchange couplings J n as- 
sumed to have spectral power density S oc l/k a . A pre- 
vious study reported some finite-size scaling evidences of 
the emergence of a phase of extended low-energy exci- 
tations [fl4| . Here, we will use a renormalization group 
technique to provide accurate estimates for the mobility- 
edge energy as a function of the degree of correlation and 
to obtain the scaling behavior governing the vanishing of 
the mobility edge in the vicinity of a = 1. Further, we 
also study the quantum diffusion of the wave-packets in 
these chains using direct integration of the motion equa- 
tions and exact diagonalization to investigate the possible 
emergence of a new dynamical regime associated with the 
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occurrence of extended low-energy spin waves for a > 1. of the Hamiltonian are therefore a linear combination of 

<fi n , i.e., $ = ^2 c n<fin, the coefficients c n satisfying the 
equation below (jlj g |l| 



II. THE MODEL AND 
RENORMALIZATION-GROUP CALCULATION 

We consider a Hamiltonian model describing a 
spin— 1/2 quantum ferromagnetic Heisenberg chain of N 
sites with random exchange couplings J„: 



H = 



N/2 
n=-N/2 



(1) 



where S n represents the quantum spin operator at site 
n and open boundary conditions are used. We take the 
exchange couplings J n connecting sites n and n + 1 to be 
correlated in such a sequencyto describe the trace of a 
fractional Brownian motion |15|, [jlj : 



N/2 
k=l 



1 1/2 



/ 2imk 
cos I — — + (fi k 



(2) 



where fc is the wave- vector of the modulations on the ran- 
dom coupling landscape and (fik are N/2 random phases 
uniformly distributed in the interval [0,27r]. The ex- 
ponent a is directly related to the Hurst exponent H 
(a = 2H + 1) of the rescaled range analysis, which de- 
scribes the self-similar character of the series and the per- 
sistent character of its increments: for a > 2 (H > 1/2) 
the increments are persistent, while for a < 2 (H < 1/2) 
they are anti-persistent. In the case of a = 2 (H = 1/2) 
the sequency of exchange couplings resembles the trace of 
the usual Brownian motion, while for a = (H = —1/2) 
one recovers the uncorrelated random exchange Heisen- 
berg model. The coupling distribution is Gaussian for 
a = but assumes a non-Gaussian form for finite a once 
the presence of long-range correlations implies in the lack 
of self-averaging and the breakdown of the central limit 
theorem. In order to avoid a vanishing exchange coupling 
we shift all couplings generated by Eq. (2) such to have 
average value < J n >= —4.5 and variance < AJ„ >= 1. 
Note that, in such case, all moments of the resulting 
distribution are finite. In order to keep the variance 
size independent, the normalization factor scales with 
the chain size. A detailed finite-size scaling analysis has 
shown how such normalization procedure is reflected in 
the main character of the one-magnon excitations p4| . 
Indeed, without such rescaling of the potential the dis- 
order width diverges for any a ^ and all states are 
expected to remain localized |ll| . 

The ground state of the system contains all spins point- 
ing in the same direction. If a spin deviation occurs at a 
site n, this excited state is described by 



S+|0>, 



(3) 



where the operator creates a spin deviation at site 
n and |0 > denotes the ground state. The eigenstates 



(J„_i + J n )c n — J„_ic„_i — J n c n+ i — 2Ec n , (4) 

where E is the excitation energy. 

In order to study the properties of one- 
magnon states, we apply a decimation 
renomalization-group technique. The method is 
based on the particular form assumed by the equation of 
motion satisfied by the Green's operator matrix elements 
[G(^)] w = (i|l/(^-H)|i)@Q: 



(E- 



mm 



S/j.,0 + J n + l i,n+ l i-l[ ( ^{ E )\n+IJ,--L,3 + J n ^ 



d n+l,n ~ u n/ 



where e£ = (J n -i + J n )/2 and J^ n+1 
After eliminating the matrix elements associated with a 
given site, the remaining set of equations of motion can 
be expressed in the same form as the original one, but 
with renormalized parameters: 
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Jn-i,n,(6) 



E-e^iE) 



-Jn- 



1,N, 



(7) 



where, after N — 1 decimations, e N denotes the renor- 
malized diagonal element at site N and Jqj/ indicates 
the effective renormalized exchange coupling connecting 
the sites and N. 

To investigate the localized/delocalized nature of the 
spin-wave modes, we compute the inverse of the excita- 
tion width or the Lyapunov coefficient 'y(E) (inverse lo- 
calization length). The Lyapunov coefficient is asymptot- 
ically related to the effective exchange coupling by M, [l9| : 



l(E) 



lim 



1 



[-ln\4% f \E)\ 



(8) 



After a linear regression of In | Jq'a/^ (E) \ versus TV we 
have a direct extrapolation of the Lyapunov coefficient 
in the thermodynamical limit. We computed j(E) for 
distinct values of the exponent a and TV = 10 5 sites. In 
addition, the density of states (DOS) was calculated by 
using the numerical Dean's method pp| |. In Fig. 1 we 
show the normalized DOS for chains with N = 10 5 sites. 
Notice that it becomes less rough as a is increased and 
its singularity at the bottom of the band is not affected 
by the imposed long-range correlation in the coupling 
constants (see insets). For a = 1.5 it consists of a non- 
fuctuating part near the bottom of the band with the 
same form as that of the pure chain (J„ = constant). 
Previous studies have pointed out that the smoothing 
of the DOS is usu ally connected with the emergence of 
delocalized states pl[ . 

In Fig. 2 we show the plot of 7 versus E for a = 
(uncorrelated random exchanges). In this case, the Lya- 
punov coefficient is finite for all energies, except for E = 
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as usual. This behavior remains qualitatively unaltered 
for < a < 1, therefore implying in the absence of ex- 
tended spin waves in this regime. For a = 0, we have 
also observed that, near the bottom of the band, the 
Lyapunov coefficient vanishes as 7 a E v with v = 1.0, 
in agreement with Ref. (h]] for probability distribution 
functions with finite < 1/J > and < 1/J 2 > (see in- 
set of Fig. 2). The picture is qualitatively different for 
a > 1. In Fig. 3(a) we plot 7 versus E for a = 1.5. The 
Lyapunov coefficient vanishes within a finite range of en- 
ergy values, thus confirming the presence of low-energy 
extended spin waves. In all chains studied with sizes 
ranging from 10 5 up to 10 6 sites, the ~f{E) curves ap- 
pear to be the same, indicating that the extended phase 
of magnons is stable in the thermodynamical limit. The 
phase diagram in the (E c , a) plane is shown in Fig. 3(b), 
with E c (given in units of A J) denoting the mobility edge 
and statistical errors are smaller than the symbol sizes. 
The data analysis (see inset) suggests that, at the vicin- 
ity of a = 1, the mobility edge depends on the correlation 
exponent as E c oc (a — l) 7 , with 7 = 1/3. 



III. SPIN- WAVE DYNAMICS 

In order to investigate the spin-wave dynamics, we 
compute the time dependence of the mean-square dis- 
placement of the wave-packet. Let us consider an excita- 
tion initially localized at site no, represented at t = by 
its eigenfunction 4> n (t = 0) = Its time evolution is 

described by the Schrodinger equation (h = 1): 



, #nft) 
1 dt 



H<l> n {t), 



(9) 



whose time-dependent wave-function can be written in 
terms of the computed eigenvectors V^' and eigenvalues 
Ej of H as @ 



Fig. 4 we show the mean squared displacement versus 
time for a = as obtained by integrating the equations 
of motion. In this case, for longer times the wave-packet 
presents a super-diffusive spread [c{t) 2 oc £ 3 / 2 ] in agree- 
ment with previous studies [Tl} | for uncorrelated random 
exchange distribution with < 1/J > finite. In Fig. 5 
we plot the data obtained by integrating the equations 
of motion for a = 0.5 and 20000 sites. As indicated, 
the initial ballistic motion extends over longer times, al- 
though a super-diffusive motion still takes place after this 
initial transient. Finally, Fig. 6 shows that for a = 1.5 
the wave-packet displays only a ballistic spread. In this 
case our calculation was performed by using both numer- 
ical integration of the equations of motion for N = 20000 
sites and exact diagonalization for N = 2000, 4000 and 
8000 sites, in which case the end effect is present (see 
inset). 

In contrast with the case of uncorrelated random ex- 
change couplings [ [j"l| , in which an asymptotic super- 
diffusive behavior sets up whenever < 1/J > is finite, we 
find a crossover with increasing a from super-diffusive to 
ballistic asymptotic regimes induced by long-range cor- 
relations in the exchange couplings. The ballistic regime 
for a > 1 can be understood following arguments sim- 
ilar to those used in Ref. pi 11 . Exploring the exact 
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mapping of the magnon problem onto a paired electronic 
one, the diffusion coefficient D can be estimated by in- 
tegrating v(k)X(k) over the extended states that effec- 
tively participate in the transport, where v(k) and X(k) 
are respectively the velocity and mean free-path of the 
excitation mode with wavenumber k. In a finite chain 
all extended modes have X(k) ~ N and travel with finite 
velocity since in the electronic problem the DOS is non- 
singular near the band center. Once there is a finite frac- 
tion of states that are dclocalizcd, the integration runs 
over a finite wavenumber and, interchanging N and t, 
the diffusion coefficient results D oc t. Consequently, the 
mean square displacement (a 2 (t)) = Dt oc t 2 confirming 
the ballistic nature of the wave-packet spread found in 
our numerical analysis. 



The same time-dependent wave-function can be obtained 
by integrating the equations of motion. The second mo- 
ment of the corresponding spatial probability distribu- 
tion is then given by 



'(*) 



(n-no) 2 cf> n (t)<l>*(t). 



(11) 



From the mean-square displacement <r 2 (t) we can es- 
timate the wave-packet spread in space at a time t. For 
any a > 0, we find ballistic behavior [a(t) 2 oc t 2 ] for ini- 
tial times, indicating that the disorder have not yet been 
realized by the spin waves p[ . In the case of uncorre- 
lated random exchange (a = 0) the self-expanded chain 
was used to minimize end effects. When the probability 
of finding the particle at the ends of the chain exceeded 
10~ 100 we added new sites thus expanding the chain. In 



IV. CONCLUSIONS 

In summary, we have studied the one-dimensional 
quantum Heisenberg ferromagnet with exchange cou- 
plings exhibiting long-range correlated disorder with 
spectral power density proportional to l/k a . By using 
a decimation renormalization-group technique we have 
found further evidences suggesting that this magnetic 
system displays a phase of extended spin waves in the 
low-energy region for a > 1 (H > 0). The mobility edge 
separating low-energy extended and high-energy local- 
ized states was shown to depend on the degree of cor- 
relation in a very special manner. Finally, through in- 
tegration of the equations of motion and exact diago- 
nalization, we have also computed the mean-square dis- 
placement of the spin- wave packet. For < a < 1, we 
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have found long-time super-diffusion, in agreement with 
previous works for uncorrelated random exchange distri- 
bution with < 1/J > finite. However, for strong corre- 
lations (a > 1) a long-time ballistic regime was numer- 
ically observed which is associated with the emergence 
of extended excitations. We believe that the reported 
results might be useful to stimulate further theoretical 



and experimental investigations of spin-wave dynamics 
on correlated ferromagnetic chains and non-periodic fer- 
romagnetic super-lattices. 
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Figure Captions 

Fig. 1 The spin-wave density of states (DOS) for chains 
with N = 10 5 sites using Dean's method. The DOS 



becomes less rough as a is increased and its singularity 
at the bottom of the band is not affected by the imposed 
long-range correlation in the coupling constants (insets). 
For a = 1.5 it consists of a non-fuctuating part near the 
bottom of the band with the same form as that of the 
pure chain (J„ = constant). 



Fig. 2 Lyapunov coefficient 7 versus energy E for a = 
(uncorrelated random exchange model) and A" = 10 5 
sites from the renormalization procedure. The Lyapunov 
coefficient is finite for non-zero energies (localized states), 
and vanishes as 7 cx E v ', E — ► 0, with v = 1.0 (inset). 



Fig. 3 (a) Lyapunov coefficient 7 versus energy E for 
a = 1.5 and N = 10 5 sites from the renormalization pro- 
cedure. The Lyapunov coefficient vanishes within a finite 
range of energy values revealing the presence of extended 
low-energy spin waves, (b) (E c ,a) phase diagram, where 
E c is the mobility edge (in units of A J) for N = 10 5 sites. 
The phase of extended spin-waves emerges for a > 1 and 
E < E c . 



Fig. 4 Mean squared displacement a 2 versus time t for 
a = from the integration of the equations of motion. 
The self-expanded chain was used to minimize end ef- 
fects. The spread of the wavepacket depicts a crossover 
from a initial ballistic (a 2 cx t 2 ) to a super-diffusive 
(er 2 oc i 3 / 2 ) behavior. 

Fig. 5 Mean squared displacement a 2 versus time t for 
a = 0.5 and N = 20000 sites from the integration of the 
equations of motion . A longer living ballistic motion 
(er 2 cx t 2 ) is found but still followed by a crossover to the 
super-diffusive regime (a 2 oc t 3 / 2 ). 

Fig. 6 Mean squared displacement a 2 versus time t for 
a = 1.5 for N = 20000 sites from the integration of 
the equations of motion. Inset: results using exact di- 
agonalization for N — 2000 (dotted line), 4000 (dashed 
line) and 8000 (dot-dashed line) sites. Ballistic behavior 
(a 2 oct 2 ) is found for all times. 



